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Abstract We discuss the possibility that the detection of gravitational waves 
emitted by compact stars may allow to constrain the MIT bag model of quark 
matter equation of state. Our results show that the combined knowledge of 
the frequency of the emitted gravitational wave and of the mass, or the radi- 
ation radius, of the source allows one to discriminate between strange stars 
and neutron stars and set stringent bounds on the bag constants. 
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What is the absolute ground state of matter? A definite answer to this 
fundamental question is still lacking. Following Bodmer's seminal paper pQ, in 
1980s Witten [5] suggested that matter consisting of degenerate up, down and 
strange quarks, called "strange" quark matter, might be bound and stable 
at zero temperature and pressure. This hypotesis, while not being directly 
verifiable in terrestrial laboratories, may be confirmed by the observation of 
"strange stars" , i.e. compact astrophysical objects entirely made of strange 
matter except, possibly, for a thin outer crust. 

Theoretical studies of structural properties of quark stars, pioneered by 
Itho over three decades ago [3], have been booming after ROSAT discovered 
the isolated pulsar RXJ1856. 5-3754 4 . The data reported in Ref. [5] seemed 
in fact to indicate that its radiation radius should be in the range 3.8 — 8.2 km, 
so that the corresponding stellar radius would be far too small compared to 
that typical of neutron stars. The results of this analysis triggered a number 
of speculations on the nature of RXJ1856. 5-3754, and it was suggested that 
it may be a strange star, though this hypotesis has now been ruled out [6]. 

Whether quark stars do exist in nature is an open question, and the 
astrophysical scenario in which they may form is still poorly understood. For 
instance, it has been argued that, besides the standard gravitational collapse, 
strange stars may form in low-mass X-ray binaries when, due to accretion, 
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matter in a neutron star core reaches sufficiently high densities to undergo a 
deconfinement phase transition to quark matter [7]. 

In this letter we discuss the possibility that detection of a gravitational 
signal emitted by a compact star, oscillating in its fundamental mode with 
frequency Vf and damping time Tf, will allow to infer whether the source is 
a neutron star or a strange star and to constrain theoretical models of the 
quark matter equation of state (EOS). 

Due to the complexity of the fundamental theory of strong interactions 
between quarks (Quantum Chromo-Dynamics, or QCD), theoretical studies 
of strange stars are necessarily based on models. The most used is the MIT 
bag model [S], in which the two main elements of QCD, namely color con- 
finement and asymptotic freedom, are implemented through the assumptions 
that: i) quarks occur in color neutral clusters confined to a finite region of 
space (the bag), the volume of which is limited by the pressure of the QCD 
vacuum (the bag constant B), and ii) residual interactions between quarks 
are weak, and can be treated in low order perturbation theory in the color 
coupling constant a s . 

The bag model parameters are the bag constant B, the quark masses m/ 
(the index / = u,d, s labels the three active quark flavors; at the densities 
relevant to our work heavier quarks do not play a role) and the running 
coupling constant a s , whose value at the energy scale /j, can be obtained 
from the renormalization group relation 

asM = (33-2N f )]n(f/A* iCD ) ■ (1) 

In the above equation Nj — 3 denotes the number of active flavors, while 
Aqcd is the QCD scale parameter, whose value is constrained to the range 
100 — 250 MeV by high energy data (see, e.g., [9]). At the scale typical of the 
quark chemical potentials Eq. (p} yields 

a. G [0.4,0.6] . (2) 

As quarks are not observable as individual particle, their masses are not 
directly measurable. However, they can be inferred from hadron properties, 
supplemented by theoretical calculations. According to the 2002 Edition of 
the Review of Particle Physics fTU], the masses of up and down quarks do 
not exceed few MeV, and can therefore be safely neglected, while the mass 
of the strange quark is much larger, its value being in the range 

m s e [80, 155] MeV . (3) 

The bag constant is subject to a much larger uncertainty. In early applications 
of the MIT bag model B, a s and m s were adjusted to fit the measured 
properties of light hadrons (spectra, magnetic moments and charge radii). 
This procedure leads to values of B that differ from one another by as much 
as a factor of ~ 6, ranging from 57.5 MeV/fm 3 [11] to 351.7 MeV/fm 3 [12], 
while the corresponding values of a s turn out to be close to or even larger 
than unity. Moreover, the strange quark masses resulting from these analyses 
are typically much larger than the upper limit given by Eq. ([3]). Using the 
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parameters determined from fits to light hadron spectra to describe bulk 
quark matter is questionable, as these spectra are known to be strongly 
affected by the details of the bag wave- functions, as well as by spurious 
contributions arising from the center of mass motion. 

The requirement that strange quark matter be absolutely stable at zero 
temperature and pressure implies that B cannot exceed the maximum value 
Bmax ~ 95MeV/fm 3 [13]. For values of B exceeding B max , a star entirely 
made of deconfmed quarks is not stable. Under these conditions quark matter 
can only occupy a fraction of the available volume and the star is said to be 
hybrid [H] . The results of Ref. [T5lfT6] , based on a state-of-the art description 
of the low-density hadronic phase, suggest that, assuming that the transition 
to quark matter proceeds through the formation of a mixed phase, hybrid 
stars can only exist for values of the bag constant up to ~ 200 MeV/fm 3 
and contain a rather small amount of quark matter. Gravitational emission 
from the hybrid star models of Ref. [15], corresponding to m s = 150 MeV, 
a s = 0.5 and B = 120 and 200 MeV/fm 3 , referred to as APRB120 and 
APRB200, respectively, has been analyzed in Ref. |17j . 

In this paper we focus on bare strange stars and consider values of the 
bag constant in the range 

B e [57, 95] MeV/fm 3 . (4) 

The problem we shall investigate is the following. Suppose that a gravita- 
tional signal is detected, which is emitted by a compact object oscillating in 
its fundamental mode (f-mode), i.e. the mode which is known to be the most 
efficient as far as gravitational radiation is concerned [T5]. We do not know 
whether the source is a neutron star or a strange star. Since the damping 
time of the f-mode , r/, is known to be of the order of a fraction of second, 
the mode excitation would correspond to a sharp peak in the energy spec- 
trum of the detected signal, emitted at the mode frequency vj which could 
be identifiable by a suitable data analysis technique. The questions we want 
to address are: 

— Does the knowledge of Vf (and/or of r/) allow one to say anything about 
the nature of the source? 

— Assuming that we can establish the star is a strange star, would these data 
allow one to set constraints on the parameters of the MIT bag model? 

To answer these questions, we compute frequency and damping time of the 
fundamental mode of strange stars, letting the parameters of the bag model 
vary in the range indicated by Eqs. {3]) and f4|, which covers the param- 
eter space allowed for bare strange stars. We consider masses in the range 
[0.7 Mq, M max ], where M max is the maximum mass allowed by each choice of 
the model parameters. We consider bare stars without a crust, as the presence 
of a crust does not affect the fundamental mode frequency in a significant 
way. 

We compare these frequencies with those computed in |17] for neutron 
stars and for hybrid stars. In [T7] NS were modeled using a set of modern EOS 
that describe matter at supranuclear densites; they are obtained within non 
relativistic nuclear many-body theory and relativistic mean field theory, that 
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model hadronic interactions in different ways, leading to different composition 
and dynamics. The hybrid stars were modeled using the EOS APR120 and 
APR200 of Ref. [15] , describing hybrid stars with a rather small admixture 
of quark matter. 
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Fig. 1 The frequency of the fundamental mode is plotted as a function of the mass 
of the star, for neutron/hybrid stars described by the EOS employed in |17| and 
indicated with the same labels (APR2, APRB200, APRB120, BBS1, G240) and 
for strange stars. The shaded region covers the range of parameters of the MIT 
bag model considered in this paper, i.e. a s £ [0.4,0.6], m a £ [80,155] MeV and 
B 6 [57,95] MeV/fm 3 . 
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Fig. 2 The damping time of the fundamental mode is plotted as a function of 
the mass of the star, as in Fig. 1, for neutron/hybrid stars and for strange stars 
(shaded region). 

The results of this comparison are shown in Figs. 1 and 2, where we plot 
Vf and the corresponding damping time r/, respectively, versus the mass of 
the star. In both figures the shaded region refers to all values of and Tt 
allowed for strange stars, assuming that the paramaters of the bag model 
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Fig. 3 The fundamental mode frequency Uf is plotted versus the gravitational 
mass, M, for different values of the bag constant and a 3 and m s varying in the 
range indicated by Eqs. ((2| and Q. 
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Fig. 4 vf is plotted versus the radiation radius, 7?oo, for different values of the 
bag constant and a 3 and m B varying in the range indicated by Eqs. ([2]) and ([3}- 
The continuous lines refer to the values of Vj for the neutron/hybrid star models 
considered in [17] . 



vary in the selected range. Continuous lines refer to the values of Vf and r/ 
computed for neutron and hybrid stars in [17j . Figure 1 shows that: 

— Strange stars cannot emit gravitational waves with Vf < 1.7 kHz for any 
values of the mass in the range we consider. 

— For masses lower than 1.8 Mq, above which no stable bare strange star can 
exist, there is a small range of frequency where neutron/hybrid stars are 
indistinguishable from strange stars. However, there is a large frequency 
region where only strange stars can emit. For instance if M — 1.2 M , 
a signal with Vj > 1.9 kHz would belong to a strange star. Note that 
the fundamental mode frequency and damping time of the hybrid stars 
considered in [IS] and [T7] and shown in Figs. 1 and 2 (EOS APR120 and 
APR200), are basically indistinguishable from that of the neutron star 
with the same low density EOS, except when the mass is close to the 
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maximum mass. This is due to the small amount of quark matter these 
hybrid stars contain. 
— Even if we do not know the mass of the star (as it is often the case for 
isolated pulsars) the knowledge of vt allows to gain information about 
the source nature; indeed, if Vf > 2.2 kHz we can reasonably exclude that 
the signal is emitted by a neutron star. 

Figure 2 contains a complementary information: for strange stars r/ is in 
general smaller than for neutron/hybrid stars. 

Thus, the next question is: assuming that we know the signal has been 
emitted by a strange star, can we constrain in some way the parameters of 
the MIT bag model ? In Fig. 3 we show to what extent this is possible. 
We plot the values of Vf allowed for strange stars versus the stellar mass, 
indicating with the same symbol the points that belong to the same value 
of the bag constant B. From this picture we see that for a given mass the 
mode frequency increases with B, and that, knowing M and Vf we would be 
able to set constraints on B much more stringent than those provided by the 
available experimental data. 

Similar information can be derived by the simultaneous knowledge of Vf 
and the radiation radius 

Roo = ; = (5) 
y/1 - 2M/R K ' 

In Fig. 4 we plot Vf as in Fig. 3, but versus Roo- In addition we plot as con- 
tinuous lines the values of Vf for the neutron/hybrid star models considered 
in [17]. The figure shows that radiation radii smaller than ~ 13 km should 
be attributed to strange stars, whereas if 13 < Roo < 15.5 Km the star can 
either be a strange or a neutron/hybrid star. Higher values of Rao can only 
belong to neutron/hybrid stars. Figure 4 further shows that the knowledge 
of Vf constrains the value of B. 

The problem whether quark stars can be discriminated from neutron stars 
using gravitational waves, which we discuss in this letter, has already been 
addressed in [19], [20] and [21]. 

In [19] Vf and r/ have been computed for strange star models obtained 
within the MIT bag model for two values of B, i.e. B = 56 and 67 MeV/fm 3 , 
m s — 150 MeV and assuming a s both vanishing and different from zero. 

In [20] strange stars have been modeled putting a s and m s to zero and 
choosing B = 75 and 137 MeV/fm 3 . 

Finally, in [5T] the values of Vf and r/ have been plotted versus the ra- 
diation radius for the following values of the parameters: a s = 0.6, m s = 
0,150,300 MeV, and B — 57 and 209 MeV/fm 3 . These values of the pa- 
rameters and the star central density were chosen to fit Roo within the range 
3.8 — 8.2 Km. However, this choice implies very small values of the stellar mass 
(ranging from ~ 0.05 to ~ 0.5 Mq) which are hard to explain within current 
evolutionary scenarios for neutron stars formation. As a consequence, they 
obtain values of Vf larger than ~9 kHz, much higher than those we consider. 

The study we propose in this paper differs from the preceeding literature 
in what it explores the entire range of allowed parameters in a systematic 
way. Our results show that the detection of a signal emitted by a compact 
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star pulsating in its fundamental mode, combined with a complementary 
information on the stellar mass or the radiation radius, would allow one to 
discriminate between neutron/hybrid stars and strange stars; in addition, 
we show that it would also be possible to constrain the bag constant to a 
range much smaller than that provided by the available data from terrestrial 
experiments. 

The fundamental mode can be excited in a variety of astrophysical pro- 
cesses, like in a glitch, in a close interaction with a companion, or after birth 
in a gravitational collapse. Recent simulations of gravitational collapse show 
that a significant fraction of the total energy emitted in gravitational waves, 
of the order of 10~ 9 — 10~ 8 M Q c 2 , is indeed emitted at the frequency of the 
f-mode [32], [53]. However, this energy is too low to be detectable by current 
interferometric antennas like VIRGO or LIGO, unless the collapse occurs in 
our galaxy; but we know that, unfortunately, the rate of collapse per galaxy 
is only of a few per hundred years. In order to detect signals emitted by more 
distant sources, we would need very sensitive, high frequency detectors, like 
EURO or EURO-XYLO, which have been considered in a preliminary as- 
sessment study some years ago (http://www.astro.cf.ac.uk/geo/euro/!. As 
discussed in [24] , this kind of detectors would be able to see signals emitted 
by oscillating stars up to the distance of the VIRGO cluster if the energy 
stored in the mode is of the order of 1CU 7 — 1CP 8 M Q c 2 , which is not too 
far from present estimates. Indeed, current numerical simulations assume ax- 
isymmetric collapse, but if asymmetries are present the emitted energy may 
be larger. And moreover rotation, which is certainly present in stars, has the 
effect of lowering the mode frequencies, enhancing detection chances. 

Thus, we can conclude that asteroseismology will become a branch of 
gravitational wave research when high frequency detectors will be operating, 
and we hope that the EURO project will be reconsidered in a not too far 
future. 

It has also to be mentioned that, although somewhat more refined dy- 
namical models have been proposed, the MIT bag model appears to provide 
a quite reasonable description of quark matter. The results discussed in Rcf. 
[2"5] show that the EOS obtained from the Nambu Jona-Lasinio (NJL) model, 
in which quark masses are dynamically generated through the appearance of 
a condensate associated with chiral symmetry breaking, is in fact similar to 
the one obtained from the bag model. Using more sophisticated models, such 
as the NJL model, may turn out to be required to describe the possible occur- 
rence of the pairing instability induced by the attractive one-gluon exchange 
interaction between quarks, leading to a color superconducting phase |25[ 
[26] . However, the relative stability of the different superconducting phases 
discussed in the literature is not yet firmly established, and their occurrence 
is expected to affect mostly transport properties and cooling, rather then the 
stellar structure, and consequently the f-mode frequency, on which we focus 
in this paper. 



References 



1. A.R. Bodmer, Phys. Rev. D 4, (1971) 1601. 



8 



2. E. Witten, Phys. Rev. D 30, (1984) 272. 

3. N. Itoh, Prog. Teor. Phys. 44, (1970) 291. 

4. F.M. Walter, S.J. Wolk and R. Neuhauser, Nature 379, (1996) 233. 

5. J.J. Drake et al. Ap. J. 572, (2002) 996. 

6. R. Turolla, S. Zane, J.J. Drake, Ap. J. 603, (2004) 265; J.E. Triimper, V. 
Burwitz, F. Haberl and V.E. Zavlin, Nucl. Phys. B (Proc. Suppl.) 132, (2004) 
560; J.A. Pons, F.M. Walter, J.M. Lattimer, M. Prakash, R. Neuhauser and 
P.H. An, Ap. J. 564, (2002) 981. 

7. K.S. Cheng and Z. G. Dai, Phys. Rev. Lett. 77, (1996) 1210. 

8. A. Chodos, R.L. Jaffe, K. Johnson, C.B. Thorne and V.F. Weiskopf, Phys. Rev. 
D 9, (1974) 3471. 

9. R.K. Ellis, W.J. Stirling, and B.R. Webber, QCD and Collider Physics (Cam- 
bridge University Press, Cambridge, 1996). 

10. K. Hagiwara, et al, Particle Data Group, Phys. Rev. D 66, (2002) 1. 

11. T. De Grand, R.L. Jaffe, K. Johnsson and J. Kiskis, Phys. Rev. D 12, (1975) 
2060. 

12. C.E. Carlson, T.H. Hansson and C. Peterson, Phys. Rev. D 27, (1983) 1556. 

13. E. Farhi and R.L. Jaffe, Phys. Rev. D 30, (1984) 2379. 

14. N.K. Glendenning, Compact Stars (Springer-Verlag, Berlin, 1997). 

15. A. Akmal, V.R. Pandharipande and D.G. Ravenhall, Phys. Rev. C 58 (1998) 
1804 

16. O. Benhar and R. Rubino, A & A 434, (2005) 247. 

17. O. Benhar, V. Ferrari and L. Gualtieri, Phys. Rev. D 70, (2004) 124015. 

18. G. Allen, N. Andersson, K.D. Kokkotas and B.F. Schutz, Phys. Rev. D 58, 
(1998) 124012. 

19. C.W. Yip, M.-C.Chu, P.T. Leung, Astroph. J. 513, (1999) 849. 

20. Y. Kojima and K. Sakata, Prog. Teor. Phys. 108, (2002) 801. 

21. H. Sotani and T. Harada, Phys. Rev. D 68, (2003) 024019; H. Sotani, K. Khori 
and T. Harada, Phys. Rev. D 69, (2004) 084008. 

22. M. Shibata, Y. Sekiguchi, Phys. Rev. D 71 (2005) 024014. 

23. H. Dimmelmeier, J. A. Font, E. Midler, A&A 393, (2002) 523. 

24. V. Ferrari, G. Miniutti, J. A. Pons, Class. Quant. Grav. 20 n. 17, (2003) S841. 

25. M. Buballa, Phys. Rep. 407 (2005) 205. 

26. M. Alford and S. Reddy, Phys. Rev. D 67, (2003) 074024. 



